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Reflection dominated X-ray spectra of Narrow Line Seyfert 
1 Galaxies: Mrk 478 and EXO 1346.2+2645 
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ABSTRACT 

Multi-epoch XMM- Newton spectra of the two Narrow Line Seyfert 1 (NLS1), Mrk 
478 and EXO 1346.2+2645, are presented. The data were fitted with different models, 
including the relativistically-blurred reflection model REFLIONX, which was found to 
give a good description of both spectra in all epochs. The two sources are reflection 
dominated with the illuminating continuum hidden from view. This can be explained 
either in terms of a corrugated disc or strong gravitational light bending effects. The 
emission from these two sources comes from a very small region (few gravitational 
radii), very close to a rapidly spinning black hole. Spectral variability analysis show 
that both sources have constant RMS spectra, and constant hardness ratios in all 
epochs, in agreement with a single spectral component. The constancy the spectrum 
between high and low states rules out alternative models such as partial covering. 

Key words: galaxies: active - galaxies:Seyfert - galaxies: nuclei - galaxies: individual: 
Mrk 478, EXO 1346.2+2645 - X-rays: galaxies. 



1 INTRODUCTION 

The first study of Narrow Line Seyf ert 1 galaxies (here- 
after NLS1) as a class was done by lOsterbrock fc Poggel 
l|l985l ). The sources were identified as having: (1) Balmer 
lines only slightly broader than the forbidden lines such as 
[OIII], [Nil] and [SII], with widths less th an 2000 km s" 1 , 
and (2) the line ratio [OIIL1A5007/H/3 < 3 l|Goodrichlll989h . 
The latter defines them as Seyfert 1 rather than 2. In addi- 
tion, they show strong Fell emission lines and iron emission 
from higher ionisation states. 

The importance of NLS1 became apparent after the dis- 
covery of the first Principle Component Analysis eigenvec- 
tor (El) of Seyfert Is. This involves anti-correlations be- 
tween the width o f the Balmer lines and the strength of 
Fell emission line dBoroson fc Green|[l992l , see a review by 
ISulentic et "alfaool ) . An additional anti-correlation was also 
found between ROSAT s pectral index (0.1- 2.4 keV) and the 
FWHM of the H/3 line (|Boller et al.lll996l h NLSls occupy 
an extreme corner of this correlation parameter space, with 
small line widths, strong Fell emission and steep X-ray pho- 
ton indices. These properties seem to be d riven by small 
black hole masses and high Eddin gton ratios |Sulentic et alj 
l200d : lBorosonll2002l : [Grupell2004r )h 

At soft X-ray energies (less than ~ 1 keV), NLS1 galax- 
ies often show an excess above the extrapolation of a simple 



E-mail: azoghbi@ast.cam.ac.uk 



power-law from higher energies. This soft excess can usually 
be fitted by a blackbody with temperatures of ~0.1-0.2 keV 
l|Gierliriski fc Don"el|2004] ). The fact that the temperature is 
constant for different black hole masses, is sugge stiv e of an 
atomic origin for the soft excess. iGierlinski fc Donel (I2004T ) 
suggested that relativistically smeared, partially ionised ab- 
sorption can lead to an apparent soft excess below ~ 0.7 
keV. The drawback of this is that in order to produce the 
smooth spectrum, extreme ou tflow velocities (~ c) are re- 
quired (|Schurch fc D one 2007). Alternatively, the soft excess 
can be explained in te rms of relativistically blurred reflection 
|Crummv et aljfeood h where the soft excess is the result of 
highly bl urred emission lines re sulting from the reflection off 
cold gas (|Ross fc Fabianll2005h . 

In the hard band (~1-10 keV), NLS1 galaxies show in 
general a steeper photon index compared to l normaV Broad 
Line Seyfert 1 , which is attributed to Compton cooling o f 
the corona that emits the hard X-rays (|Brandt et al.lll997T l. 
A correlation was found between the spectral ind ex and 
the Eddington ratio in sam ples that include NLS1 l|Grupel 
l2004l ; ls"hemmer et al.ll2008T ), and points towards the fact that 
NLS1 accrete at high rates. Another feature that is some- 
times observed in NLS1 is a spectral drop around 7 keV , first 
observed in the source 1H 0707-495 (|Boller et al.ll2002f). and 
later in other sources (e.g. IRAS 13224-3809. iBoller et aT 



2005tl . The feature appears to change in time IjGallo et al 



20041 for 1H 0707-495) and is not seen in all NLS1 AGN. Of- 
ten it may appear as spectral curvature around 7 keV (e.g 
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Table 1. XMM-Newton observation log. 



Source z Galactic absorption Obs. Obs. Date PN net 6 Average PN (0.2-10 keV) 

Name Njj(x 10 20 cm — 2 ) a num. ID exposure (ks) count rate (counts/s) 

Mrk 478 0.079 1.03 1 0107660201 23/12/2001 10.1 7.5 

2 0005010101 01/01/2001 14.5 5.2 

3 0005010201 04/01/2001 7.2 5.3 

4 0005010301 07/01/2001 18.2 2.8 
EXO 1346.2+2645 0.0589 1.19 1 0097820101 26/06/2000 41.9 1.2 

2 0109070201 13/01/2003 45.7 0.3 



a from Dickcv & Lockman (1990). 
b after background flare subtraction. 

lGalloll2006h . The explanation of this drop is still not clear. 
Two models that have been suggested, both giving good fits 
to the data, are partial covering and reflection. 

In the partial covering model, the central source is cov- 
ered by a patchy absorbing cloud. The absorption produces 
the sharp edge at high energies, while a fraction of the emis- 
sion leaks through to produce the soft part of the spectrum. 
This, however, requires an excess ively large ir on abundance 
of ~ 30 x solar (|Boller et alj|200"3 , see however iTanaka et al.l 
12004 ) . It should be noted that partial covering provides good 
fit for other so urces where the spe ctral drop is not seen or not 
as sharp (e.g. iGrupe et al.l 120071 ) . On the other hand, par- 
tially ionised reflection off an optically thick m aterial (an ac- 
cretio n disc) can explain the spectra very well (jFabian et al.l 
120041 ). In this model, the spectral drop around 7 keV is the 
blue wing of a relativistically broadened iron line. 

In this paper we explore the spectra of two NLS1 using 
multi-epoch observations with XMM-Newton. The paper is 
organised as follows. In section 2, a description of the ob- 
servations and data reduction is presented. Spectral analysis 
and results are given in section 3, with subsections discussing 
individual sources. A study of the spectral variability is pre- 
sented in section 4, and the implications are discussed in 
section 5. 



The pattern statistics were checked for possible pileup (us- 
ing epatplot in SAS), this was found to be present in the 
first Mrk478 observation. To correct for it, the central region 
of the source was excluded from the analysis. The pattern 
selection used is pattern ^ 4 for EPIC-PN and pattern 
«S 12 for EPIC-MOS. 

Sources spectra were extracted from circular regions 35 arc- 
sec across centred on the source, and the background from 
regions on the same chip and away from the source. The 
spectra were then grouped such that each bin has a mini- 
mum of 20 counts. EPIC responce matrices were generated 
using the tasks rmfgen and ARFGEN in SAS. It should be 
noted that EXO 1346.2+2645 is close in sky position to the 
cluster Abell 1795, and different background regions have 
been tested to make sure the spectrum of the source is 
not contaminated. Also, for the second observation of EXO 
1346.2+2645 (0109070201), the calibration reported the cen- 
tral pixels of the source as bad so they were excluded from 
the anal ysis. Spectral fitting was performed using xspec 
V12.4.0 (|A rnau dl ll996h . All quoted errors on the model pa- 
rameters correspond to a 90 % confidence level for one in- 
teresting parameter (i.e. a A\ 2 = 2.71) in the rest frame of 
the source. 



2 OBSERVATIONS & DATA REDUCTION 

The data presented here for the two sources Mrk 478 and 
EXO 1346.2+2645 consist of several observations which 
are detailed in Table [T] Some of the data have been pub- 
lished before as par t of samples (e.g. iPorquet et al. 2004, 
I Crummy et al. l l2006l . lGallcll2006h , but the majority is pre- 
sented properly for the first time here. 

The Observation Data Files (ODFs) were retrieved from 
the archive and processed in the standard way using XMM- 
Newton Science Analysis System (SAS V7.1.0). The EPIC 
cameras were operated in large window mode in the first Mrk 
478 observation (0107660201) and in small window mode in 
the last three. EXO 1346.2+2645 was observed in full frame 
mode. The data were corrected for hot, dead and flicker- 
ing pixels. Good Time Intervals (GTI) were produced in the 
standard way (SAS user guide □). The most affected were 
the first and third observations of Mrk 478 (see Table [TJ. 

1 http:/ /xmm.esa.int/external/xmm_user_support/document- 
ation/sas_usg/USG/node57.html 



3 SPECTRAL ANALYSIS 

As a start, a simple power-law was fitted to each source in 
the energy range 0.3-10 keV. Galactic absorption was in- 
cluded using the column densities shown in Table Q] With 
multi-epoch observations, all fit parameters were tied except 
for a multiplicative constant. The PN and MOS data were 
fitted separately and they gave consistent results. In what 
follows (unless otherwise stated) only the PN data are used. 
The simple power-law fit was not adequate, giving \ 2 /d.o.f 
3236/1760 and 903/787 for Mrk 478 and EXO 1346.2+2645 
respectively. The soft band fitted well leaving most of the 
residuals at higher energies (> 4 keV). 

The data were then fitted in the range 1.5-4.5 and 8.5-10 
keV. This excludes the soft band where the soft excess dom- 
inates the spectrum and also excludes the region around the 
iron complex. A good fit was found with xV rf - -/ 674/685 
and 222/289 for the two sources respectively. Fig. Q] shows 
the ratio of data/model when extrapolated over the 0.3-10 
keV energy range. This clearly shows the soft excess emis- 
sion below ~ 1.5 keV. Also, it shows residuals around ~ 6 — 7 
keV that might be caused by the presence of iron emission. 

To find a phenomenological description of the data, each 
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Figure 1. Data/model ratio of an absorbed power-law fit to the EPIC pn in the range 1.5-4.5, 8.5-10 keV extrapolated over the entire 
0.3 to 10 keV for the two sources (the power-law was multiplied by a constant to account for flux variations). Different colours represent 
the different observations as shown in Tabic [T] Data have been re- binned for display. 




Energy (keV) 

Figure 2. The best fit models for Mrk 478 in the range 0.3- 
10 keV with residuals in terms of sigmas. top: a single com- 
ponent reflection-dominated model, bottom: a partial covering 
model with a blackbody and a partially covered power-law. Both 
models include Galactic absorption, and the different colours rep- 
resent the four observations available. Data have been re-binned 
for display. 



source will be looked at individually starting by adding a 
multi-colour disc blackbody (diskbb in xspec). 



3.1 Mrk 478 

Adding the blackbody component improved the quality of 
the fit significantly (A^ 2 =1041 for 2 d.o.f) compared to 
the simple power-law, the temperature of the blackbody 
was 123^3 eV with a power-law index of F — 2.41 ± 0.02 
(compared to 89 ± 0.12 eV and 1.6 ± 0.05 respectively 
found with the ASC A data fitted between 0.6 and 10 keV, 
IVaughan et al.|[l999l ). Allowing the temperature and Y to 



change between observations changes their values by 5 and 
3 percent respectively, and improves the fits by A% 2 ~ 
17 per additional free parameter. This might imply some 
minor variability in the source between observations. This 
fit, however, leaves some residuals at high energies, which 
can be fitted by a Gaussian line with rest frame energy of 



E 



7.03±°;$l keV and a 



1.59 4 



-+0.29 



keV, and an equiva- 



lent width of EW 
This broad line implies that relativistic e ffects are at 
play. Relativistic ally broade ned lines diskline (|Fabian et al.l 
ll989T ) and Laor (|Laorlll99lf l emitted around a Schwarzschild 
and Kerr black hole, respectively, both, gave comparable 
fits (reduced \ 2 of 1.07). The inner radii are r- ln = 6.01q'q, 
rjn = 8.2l 2 'g gravitational radii (r g — GMbh/c 2 ), and line 
energies are 6.20^054, 6.491q ^3 keV, respectively, consis- 
tent with iron emission. These arise from the reflection of 
the incident power-law off the accretion disc, which implies 
the existence of a reflection continuum that has not been 
accounted for so far. 

To model it, the high resolution version of the constant 
density reflection model reflionx of lRoss fc Fabiar] (|2005l ) 
was used, blurred with kdblur kernel of lLaorl I 199ll) to ac- 
count for relativistic effects in the vicinity of the black hole. 
In this model, an optically thick disc is illuminated with a 
power-law, producing florescence lines and continuum emis- 
sion. Its parameters include the Fe abundance, ionisation 
parameter £ (defined as £ = 4nF to t/nH where Ftot is the 
total illuminating flux, and nn is the density of the reflec- 
tor) and the incident power-law index F, while the kdblur 
assumes a power-law emissivity and sharp inner and outer 
radii of the accretion disc0 (n n and r out ). The outer radius 
was frozen at 400 r g , as very little emission is expected out- 
side this radius, while all other parameters were allowed to 
vary. The fit was very good with reduced \ 2 °f 1-017 for 
1752 d.o.f. 

In order to produce the smooth soft excess, the emis- 
sion has to be highly blurred coming from a region very 
close to the black hole. The fits gave an inner radius of 
ri n = 1.231q'qq r g (the inner most stable circular orbit of 
a Kerr black hole), with a very narrow emission region im- 
plied by the emissivity index peaking at the maximum value 
allowed (lO.Oj 



^2) 



phabs(kdblur(atable{reflionx.mod})-|-powerlaw) in XSPEC. 



4 A. Zoghbi et. al 

Table 2. The parameters of the best-fit reflection models. 



Source 


Obs. 
Num 


Em. 
Ind. 




Incl.(°) 


Fc abun. 
( X solar) 


r" 


£ a (erg cm s 1 ) 


Rcf. Frac. 6 


x 2 /d.o.f 


Mrk 478 


1 

2 
3 
4 


6 C 


1 4 7 +0.08 
-0.02 


34 ± 2 


0.36 ± 0.03 


2.43 ± 0.02 
2.35 ± 0.02 

2 47+°'02 
ZMI -0.04 

2.29 ± 0.01 


1000±^ 8 

"68±fo 4 4 
746+^° 


0.95 ± 0.05 
0.81 ± 0.09 
0.88 ± 0.11 
0.89 ± 0.11 


1816/1754 


EXO 1346 


1 

2 




1 6 +0 - 6 


36 ± 8 


34+ 05 
0.03 


9 17 +0.02 
z,± ' -0.03 


1202117"; 


0.83 ± 0.09 
0.96 ± 0.04 


1255/1306 



a for Mrk 478, these parameters were allowed to vary between observations, and were tied for EXO 1346.2+2645 
b calculated as the ratio of the flux in the reflection component to the total flux. 
c parameter is frozen. 




N H (1(F) 



Figure 3. A contour plot for the 90% limits of the column density 
and covering fraction for the partial covering model for Mrk 478. 
The cross marks the best values for the combined data. The four 
colours represent the four epochs. 



The fact that some parameters (ri n and emissivity index) 
peak at the limits might be problematic, freezing the emis- 
sivity index at a lower value of 9 and 6 does not effect 
the other parameters, giving inner radii of 1.3llQ g r g and 
l-47lo 02 r g respectively. The top panel of Fig. [2] shows the 
data and residuals for the reflection model. The fit parame- 
ters are presented in Table [5] 

The normalisation of the power-law has a zero lower 
limit, which might indicate that the power-law is not re- 
quired. A reflection-only model fits the data equally well (re- 
duced x 2 /d.o.f of 1.017/1753). The fit parameters remained 
almost the same, with the values of F and £ changing by 
about 2 and 20 percent respectively when the power-law is 
removed (comparable to the errors in those values). Limits 
on the reflection fraction are shown in Table [2] 

The other model discussed in the literature is partial 
covering. A good fit to the data was found using a continuum 
of a power-law and a blackbody with a single partial covering 
component (r = 2.66 ± 0.02, and kT = 112 ± 20 keV). The 
reduced \ 2 was 1-006 for 1754 d.o.f. The data/residuals plot 
is shown in bottom panel of Fig. [2] All parameters were tied 
between observations except for the power-law index and the 
normalisations. The model has an absorption column den- 
sity of 2.7ii"i x 10 23 cm" 2 , hon abundance of 0.2+°; \ x solar 
and a covering factor of 0.54±0.03. Allowing the parameters 
to vary between observations (column density and covering 
fraction) only changes the size of the errors reflecting the 
quality of the individual observations. Fig. [3] shows the 90 
percent contour limits on the values of the column density 
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Figure 4. The best fit model for EXO 1346.2+2645 in the range 
0.3-10 keV with a single component reflection-dominated model, 
with residuals in units of sigma. Data have been re-binned for 
display, colours represent different observations (and detectors: 
PN and MOS) inset: The model used to fit the data. 

and covering fraction. The coverer parameters remain con- 
stant between epochs. 

3.2 EXO 1346.2+2645 

Fitting the data with a power-law and a blackbody gave a 
good fit with a reduced \ 2 °f 1-03 for 789 d.o.f (blackbody 
temperature: 157 ± 7 eV, and power-law index F = 2.23 ± 
0.04). Some residuals at high energies can be accounted for 
by adding a Gaussian line (E= 6.781q 34 keV, with equiva- 
lent width of l.l±%[l keV), or a Laor line (E= 6A±° l keV 
with equivalent width of 1.4 ± 0.6 keV, n n = 1.3jl r g , and 
emissivity index of 2.3l 3 +). Also, this high energy excess can 
be fitted with a partial covering model (reduced \ 2 = 1-0 for 
785 d.o.f), the required column density is 6.2lg'2 x 10 23 cm -2 
with a covering fraction of 0.65±g;?|, the iron abundance is 
2-3jl^ 4 6 x Solar. The high upper limit on the iron abundance 
indicates the data can not constrain it keV. Freezing the iron 
abundance at the solar value changes the column density to 
ll.Olg^ x 10 23 and the covering fraction to 0.74±°:32- For all 
fits to the high energy end, the value of F does not change 
significantly. 

Due to the fact that this source is faint, data from EPIC- 
MOS detectors were included in the analysis to improve 
the statistics. Fitting the whole spectrum with a reflec- 
tion model gave an equally good fit with parameters con- 
sistent with the Laor line model. The power-law was also 
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Figure 5. Background subtracted light curves for Mrk 478 (left) and EXO 1346 (right), produced with time bins of 200 seconds. The 
vertical lines seperate the different observations available. 




Figure 6. Fractional rms variability of Mrk 478, produced with 
time bins of 400 seconds. Black circles: high flux state (Obs ID: 
0107660201). Red squares: low state (Obs ID: 0005010301) 
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Figure 7. Hardness ratio for Mrk 478 in both high and low states 
(Top), with the corresponding light curve (Bottom). The hardness 
ratio was calculated as (H-S) / (H+S), where H and S are the count 
rates in the 0.3-2 and 3-10 keV bands respectively. 



not required in this source and the data can be fitted by a 
reflection-only model. Similar to Mrk 478, the outer radius 
was fixed at 400 r g , while all other parameters were allowed 
to vary. The fit gave an inner radius of rt n = 1.6lo'3r 9 , 
and emissivity index of 3.9jlJ' \ with an ionisation parame- 
ter of 1202^^2 er S cm s_1 - The best fitting model is shown 
in Fig. H 



4 SPECTRAL VARIABILITY 

In order to gain more insight into the nature of these two 
sources, light curves in different energy bands have been 
analysed. Fig. [S] shows the light curves in the energy range 
0.2-10 keV for the two sources. It shows that the flux of both 
sources change by a factor of ~ 5 between observations that 
are separated by at least a year, in addition to some small 
variability within single observations. 

To study the variability in different ene rgy bands, we 
used the fractional variability amplitude -F var (|Edelson et al.l 
120021 ) and hardness ra tios. FVar was c a lculat ed following the 
methods described in lEdelson et ail i|2002l ) for six energy 
bins. Fig. [6] shows a plot of F vax for Mrk 478 in the "low" 
and "high" states, taken during the first and last observa- 
tion (see light curve in Fig. [5}. 

It is clear that fractional variability is consistent with a con- 
stant value within each observation and between the two 
states. The best fitted constant lines have values of 12 and 



13 percent for the high and low state, respectively. The rms 
spectra of AGN often have a peak at energies of about ~ 1 
keV, and these ca n be fully explained by a two compo- 
nent spectrum (e.g. IVaughan fc Fab ian 2004) , with one con- 
stant reflection dominated component (RDC) and a variable 
power- law component (PLC). The fact that the rms spec- 
trum is constant is in full agreement with the spectrum con- 
sisting of only one reflection component as infered from the 
spectral fitting. 

To investigate this further, Fig. [7] shows the hardness 
ratio for the low and high states (first and last observations 
respectively), along with the corresponding count rates. The 
hardness ratio was calculated using (H-S)/(H+S), with H 
and S being the count rates in the hard (3-10 keV) and 
soft (0.3-2 keV) bands respectively. It remains constant and 
does not respond to the flux changes. This again supports 
the one-component model found earlier. 

A similar trend was found in the other source (EXO 
1346), the rms spectrum is shown in figure [H] The fractional 
variability is constant in both high and low states, with the 
source being slightly more variable in the high state (black 
circles in Fig. [5} . 



5 DISCUSSION 

It is clear that a simple power-law fit leaves residuals both 
at soft and hard energies. The soft excess can be fitted with 



6 A. Zoghbi et. al 




Figure 8. Fractional rms variability of EXO 1345.2+2645, pro- 
duced with time bins of 400 seconds. Black circles: high flux state 
(Obs ID: 0097820101), and red squares are for the low state (Obs 
ID: 0109070201). 



clear that this feature is not seen in the sources presented 
here, and both models give an equally good fit to the data, 
but based on the spectral variability, the reflection model is 
favoured. The sources show flux variability of up to a factor 
of 5 between epochs, with no change in the spectrum. In a 
source covered by a line-of-sight clouds, this requires that 
the observed variability is due to variations in the intrinsic 
source not the clouds covering it (because the parameters 
describing the coverer did not change). By causality argu- 
ment, the intrinsic source is less than a lighthour across, 
and the partial coverer needs to be even smaller, yet has 
to give the same covering fraction years later, this seems 
completely implausible to us. The reflection model on the 
other hand, fits the whole spectrum vey well and is physi- 
cally motivated, with the changes between epochs explained 
by changes in the ionisation of the reflector. 



a blackbody, with temperatures of 123 jl§Q and I57I34 eV 
for the two sources in this study, consistent with other stud- 
ies that showed that the observed range of temperatures in 
these sources is limited to 0.1-0.2 keV (|Gierliriski fc Done] 
|2004 ICrummv et all I2006D . these energies are high com- 
pared to wh at is expected from a stan dard thin accretion 
disc model |Shakura fc Svunvaevl [l973l ) . A Laor line pro- 
vides a good fit for the high energy residuals. However it is 
clear that a reflection model can fit the whole X-ray range 
in the these sources. 

In terms of fitting parameters, in order to produce the 
smooth shape of the soft excess, the spectrum has to be 
highly blurred, implying emission from very close to the 
black hole. The fits presented here all point towards that, 
wi th inner radii of l ess th an 2.5r g , similar to what was found 
bv lCrummv et all |2006) for a sample of 21 sources. If the 
values of r g found here are correct, then this not only implies 
that the black holes in these sources are rapidly spinning, 
but also, that strong gravitational light bending must be 
occurring. 

For both sources studied here, the spectra are domi- 
nated by reflection with no power-law required, which im- 
plies that the illuminating continuum is hidden from view. 
This is fully supported by the spectral variability analysis, 
that showed that the whole spectrum is varying in the same 
way, implying a single spectral component. 
Although a reflection-dominated spectrum implies a devia- 
tion from symmetric geometry, it c an be easily explain ed by 
a corrugated accretion flow model ()Fabian et alJ l2002). The 
disc material is clumped into deep sheets or rings caused 
by disc instabilities, surrounded by a hot X-ray emitting 
corona. The hard X-ray continuum source is hidden for most 
observers and only the reflection off the sheet is observed. Al- 
ternatively, it can be explained in terms of str ong light bend- 
ing effects in the vicinity of t he black hole ()Miniutti et all 
l2003l . lMiniutti fc Fabianf 20041 ). In this model, emission from 
the illuminating source is bent back towards the black hole, 
while the reflected photons can easily escape towards the 
observer. 

There has been lately a controversy over the exlanation 
of the spectral features of some N LS1. In particular th e sharp 
drop at ~ 7 keV in 1H0707-495 (jBoller et al.ll2002l ). which 
was explained either by reflection or partial covering. It is 
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